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Metastable quarks of 4th generation are predicted in the framework of heterotic 
string phenomenology. Their presence in heavy stable hadrons are usually strongly 
constrained; however their hidden compositions in Heavy doubly charged baryons 
here considered are found to be still allowable: we studied their primordial quark 
production in the early Universe, their freezing into cosmic Heavy hadrons, their 
later annihilation into cosmic ray as well as their relic presence in our Universe and 
among us on Earth. We discuss also their possible production in present or future 
accelerators. Indeed if the lightest quarks and antiquarks of the 4th generation are 
stored in doubly charged baryons and neutral mesons, their lifetime can exceed the 
age of the Universe; the existence of such an anomalous Helium-like (and neutral 
Pion-like) stable particles may escape present experimental limits, while being close 
to present and future experimental test. On the contrary primordial abundance of 
lightest hadrons of the 4th generation with charge +1 can not decrease below the 
experimental upper limits on anomalous hydrogen and therefore (if stable) it is 
excluded. While 4th quark hadrons are rare, their presence may play different and 
surprising role in cosmic rays, muon and neutrino fluxes and cosmic electromagnetic 
spectra. Most of these traces are tiny, but nearly detectable. 

Key words: cosmology, cosmic rays, 4th generation, anomalous helium 
PACS: 98.80.Cq, 04.70.Bw, 04.70.Dy, 98.80.Hw(?) 



1 Introduction 



The problem of existence of quarks and leptons of 4th generation is among the 
most important in the modern high energy physics. The natural extension of 
the Standard model leads in the heterotic string phenomenology to the predic- 
tion of fourth generation of quarks and leptons [1,2] with a stable massive 4th 
neutrino [3,4,5,6]. The comparison between the rank of the unifying group Eq 
(r = 6) and the rank of the Standard model (r = 4) implies the existence of 
new conserved charges and new (possibly strict) gauge symmetries. New strict 
gauge U(l) symmetry (similar to U(l) symmetry of electrodynamics) is possi- 
ble, if it is ascribed to the fermions of 4th generation. This hypothesis explains 
the difference between the three known types of neutrinos and neutrino of 4th 
generation. The latter possesses new gauge charge and, being Dirac particle, 
can not have small Majorana mass due to see saw mechanism. If the 4th neu- 
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trino is the lightest particle of the 4th quark-lepton family, strict conservation 
of the new charge makes massive 4th neutrino to be absolutely stable. Follow- 
ing this hypothesis [1] quarks and leptons of 4th generation are the source of 
new long range interaction, (which we'll call further y-interaction) similar to 
the electromagnetic interaction of ordinary charged particles. 

Recent analysis [7] of precision data on the Standard model parameters has 
taken into account possible virtual contributions of 4th generation particles. 
It was shown that 4th quark-lepton generation is not excluded if 4th neutrino, 
being Dirac and (quasi-)stable, has a mass about 50 GeV (47-50 GeV is la 
interval, 46.3-75 GeV is 2a interval) [7] and other 4th generation fermions 
satisfy their direct experimental constrains (above 80-220 GeV) . We will con- 
centrate our attention here on given values of mass for 4th generation fermions, 
assuming that the mass of the lightest quark is about 250 GeV. 

If the lightest quark Q of 4th generation possess new strictly conserved gauge 
charge, it can decay only to 4th generation leptons owing to GUT-type in- 
teractions, what makes it sufficiently long living. If its lifetime exceeds the 
age of the Universe, primordial Q-quark (and Q-quark) hadrons should be 
present in the modern matter. If this lifetime is less than the age of the Uni- 
verse, there should be no primordial 4th generation quarks, but they can be 
produced in cosmic ray interactions and be present in cosmic ray fluxes. The 
assumed masses for these quarks make their search a challenge for the present 
and future accelerators. 

In principle, the composition of the lightest baryon of 4th generation can be: 
(Uuu) (charge +2) 1 , (Uud) (charge +1), (Udd) (charge 0) or (if mo < f^u) 
(Duu) (charge +l),(Dud) (charge 0), (Ddd) (charge —1) and the correspond- 
ing lightest mesons (Uu) (charge 0), (Ud) (charge —1) or (Du) (charge +1) 
and (Dd) (charge 0). In the present paper we mainly follow the assumption 
that vtld > mjj, which excludes the baryon with the negative charge, while 
the lightest (Ud) meson seems to be excluded by the quark model arguments. 
These arguments also exclude neutral (Udd) baryon as the lightest. It leaves 
theoretically favorable (Uud) and theoretically less favorable (Uuu) as only 
candidates for lightest baryon of 4th generation with the only possibility of 
neutral (Uu) as the lightest meson. 

We analyze the mechanisms of production of metastable Q (and Q) hadrons 
in the early Universe, cosmic rays and accelerators and reveal the possible 
signatures of their existence. We'll show that though it may be more prob- 
able that the lightest hadron of the 4th generation has the charge +1 and 
Q-quark should have decay lifetime less than the age of the Universe to avoid 
the contradiction with the negative results of the anomalous hydrogen search, 



One can call this particle A^ + . 
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the theoretically less favorable case of lightest (Uuu) baryon and (Uu) meson, 
mostly considered in the present paper, deserves special interest. Such doubly 
charged and neutral primordial hadrons can be sufficiently stable to be present 
in the modern Universe and still remain elusive for their direct searches. On 
the other hand the set of direct and indirect effects of their existence provides 
the test in cosmic ray and underground experiments which can be decisive for 
this hypothesis. In view of the interesting possibility revealed for Q-quark to 
be sufficiently stable we mainly use throughout the paper the examples, corre- 
sponding to the case, when [/-quark of 4th generation is lighter than D-quark. 
However the calculations of primordial abundance, cosmic ray production and 
effects of heavy quark annihilation bear general character and are valid for 
any charge assignment for Q and its lightest hadrons. 

In the next sections we consider U and U hadrons presence since early Uni- 
verse (Sec. 2) up to their early freezing (details in App. 1-4) and later galactic 
clustering, their annihilation and survival (Sec. 3, Appendix 5), their blow- 
ing wind onto Earth (Sec. 4) and their further annihilation and relic presence 
in cosmic rays (sec. 5), polluted matter, rarest upgoing muons and possible 
Heavy quark UHECR showers influence (App. 5-8). We consider (Sec. 6) pos- 
sible Heavy quark signatures in Tevatron and LHC. In the final Section 7 we 
summarize our present study on the novel cosmic and astrophysical scenario 
on Heavy hadrons, possible hidden and stored just here among us. 



2 Primordial [/-hadrons from Big Bang Universe 

In the early Universe at temperatures highly above their masses fermions of 4th 
generation were in thermodynamics equilibrium with relativistic plasma. Strict 
conservation of new gauge U(l) charge (y-charge) implies y-charge symmetry 
of their distribution. When in the course of expansion the temperature T falls 
down below the mass of the lightest [/-quark, to, equilibrium concentration of 
quark- ant iquark pairs of 4th generation is given by 



where g± = 6 is the effective number of their spin degrees of freedom. We use 
the units H — c — k — 1 throughout this paper. 

At the freezing out temperature Tf the rate of expansion exceeds the rate of 
annihilation to gluons UU — > gg and UU pairs are frozen out. The frozen out 
concentration at T ~ Tf ~ to/30 (See Appendix 1) is given by 




(1) 



2.5 • 10 
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m 
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Even this value of primordial concentration of £/-quarks with the mass m = 
250GeV would lead to the contribution into the modern density, which is 
by an order of magnitude less, than the baryonic density, so that [/-quarks 
can not play a significant dynamical role in the modern Universe. The actual 
value of primordial [/-particle concentration should be much smaller due to 
radiative and hadronic recombination, which reduce the abundance of frozen 
out [/-particles. 

Radiation of mutually attracting y-charges of U and U can lead to formation 
of UU systems bound by y-attraction, in which U and U rapidly annihilate. 
The rate of this radiative binding is calculated in Appendix 2 and is given by 

(av) * 6 ■ 10-13(^)2(^)9/10(250^ 1/10 cm3 

x 1 v l/30 ; K T ' K m ' s K ' 



The decrease of [Z-hadron abundance owing to UU recombination is governed 
by the equation 

^ = ~r\ ■ s ■ (av) , (4) 



where s is the entropy density (see Appendix 1). 

Its solution at Tf > T (see Appendix 2) turns to be independent on the frozen 
out concentration (2) and is given by (see Appendix 2) 

r 4 « 0.013 ( ^V /10 « 4 • IP' 16 g- {ML!} 2 , (5) 

\T f j a 2 y m Pl 250GeVV % J ' K J 



where a y ~ 1/30 is the running constant of y-interaction. 

After QCD phase transition at T = Tqcd ~ 150MeV quarks of 4th generation 
combine with light quarks into [/-hadrons. In baryon asymmetrical Universe 
only excessive valence quarks should enter such hadrons, so that [/-baryons 
and [7-mesons are formed. The details of U- and £7-quark hadronization are 
discussed in Appendix 3 and following this discussion we consider further 
mainly doubly charged (JJuu) baryon and neutral (Uu) meson 2 . 

2 The total charge of (Uuu) and (Uu) is +2, what looks like a net electric charge 
asymmetry of the Universe. However the net charge of the Universe is vanishing 
due to the presence of two corresponding excessive electrons. This pair of electrons 
compensates the positive charge of three light u quarks, and the mechanism of this 
compensation is provided by the usual generation of baryon (and lepton) asymmetry 
of the Universe. 
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As it was revealed in [1,2] in the collisions of such mesons and baryons re- 
combination of U and U into unstable (UU) "charmonium -like" state can 
take place, thus successively reducing the [/-hadron abundance. Hadronic re- 
combination should take place even in the absence of long range ^/-interaction 
of [/-particles. So, we give first the result without the account of radiative 
recombination induced by this interaction. 

The uncertainties in the estimation of hadronic recombination rate are dis- 
cussed in Appendix 4. The maximal estimation for the reaction rate of recom- 
bination (av) is given by 



1 _ ^-i6 cm: 



3 



(av) T « 6 • 10" ib (6) 

m% s 



or by 



(av) ~ i- w 2 • 10~ 17 —. (7) 

m 2 p s 



The minimal realistic estimation (See Appendix 4) gives 

(av) « 0.4 • (T e// ro 3 )- 1/2 (3 + log (T QCD /T eff )), (8) 

where T e // = max{T, atyin^}. 

Solution of Eq.(4) for (av) from the Eq.(6) is given by 
Case A 

n = ^= — — - » l.o • io- 20 (9) 

1 _|_ r W 9QCD mpi Iqcd 

L "T Y 45 myr mir 



and it is (^f-) ~ 30 times larger for (av) from the Eq.(7): 
Case B 



r 4 = f lL= » 3.0 • IO" 19 . (10) 

1 _l_ t- l ^9QCD mpi iQCD 



For the minimal estimation of recombination rate (8) the solution of Eq.(4) 
has the form 

r 4 = -J2= == » 5 • 10- 16 ( ™ ) 3 / 2 , (11) 
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where in all the cases r is given by Eq.(2) and Qqcd ~ 15 (see Appendix 3). 
These solutions are independent on the actual initial value of r 4 = r , if before 
QCD phase transition it was of the order of (2). We neglect in our estimations 
possible effects of recombination in the intermediate period, when QCD phase 
transition proceeds. 

With the account for radiative recombination it is the value (5) that should 
be taken at T ~ Tqcd as r in the solutions (9), (10) and (11) for the results 
of hadronic recombination. Such account does not change the values (9) and 
(10), which are still with the high accuracy independent on r . However, for 
the minimal estimation of the recombination rate (8) the result of hadronic 
recombination is modified and reads 

Case C 

n ~ 2 • 10~ 16 ( ™ f' 2 . (12) 
V 250GeW V ; 



The existence of new massless U(l) gauge boson (^/-photon) implies the pres- 
ence of primordial thermal y-photon background in the Universe. Such back- 
ground should be in equilibrium with ordinary plasma and radiation until the 
lightest particle bearing y-charge (4th neutrino) freezes out. For the accepted 
value of 4th neutrino mass (> 50GeV) 4th neutrino freezing out and corre- 
spondingly decoupling of |/-photons takes place before the QCD phase transi- 
tion, when the total number of effective degrees of freedom is sufficiently large 
to suppress the effects of y-photon background in the period of Big Bang nu- 
cleosynthesis. This background does not interact with nucleons and does not 
influence the BBN reactions rate, while the suppression of y-photon energy 
density leads to insignificant effect in the speeding up cosmological expansion 
rate in the BBN period. In the framework of the present consideration the 
existence of primordial y-photons does not play any significant role in the 
successive evolution of [7-hadrons. 

In the period of recombination of nuclei with electrons charged [7-baryons 
recombine with electrons to form atoms of anomalous isotopes. The substan- 
tial (up to 10 orders of magnitude) excess of electron number density over 
the number density of primordial [/-baryons makes virtually all fJ-baryons to 
form atoms. The cosmological abundance of free charged fJ-baryons is to be 
exponentially small after recombination. If the lightest is (Uuu) baryon with 
electric charge +2, atoms of anomalous He are formed, [7-hadrons with charge 
+1 form atoms of anomalous hydrogen. 

These atoms, having atomic cross sections of interaction with matter, partici- 
pate then in formation of astrophysical bodies, when galaxies are formed. One 
might assume that neutral ?7-hadrons, having nuclear interactions and being 
clustered in galaxies, should not follow matter in formation of stars and plan- 
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ets, and that one can expect suppression of their concentration in such bodies. 
However the existence of Coulomb-like y-attraction will make them to obey 
the condition of neutrality in respect to the y-charge. Therefore owing to neu- 
trality condition the number densities of U- and £7-hadrons in astrophysical 
bodies should be equal 3 . 



3 Evolution of [7-hadron content in galactic matter 



In the astrophysical body with atomic number density n the initial ?7-hadron 
abundance riuo — r u • n can decrease with time due to UU recombination. 
Under the neutrality condition 

n v = rifj 

the relative ?7-hadron abundance r = njj/n = n^/n is governed by the equa- 
tion 

— = — r 2 ■ n ■ (av) . (13) 
at 



The solution of this equation is given by 

r = Tu - - . (14) 

1 + r u ■ n ■ (av) -t 



If 



n 



{av) ■ t > 



(15) 



the solution (14) takes the form 



n • (av) • t 



3 In this sense there is a very remarkable and obvious difference between common 
(light-quark) baryon matter and the 4th quark matter: there is no room for any 
lepton-baryon asymmetry as in our matter. The same asymmetry leads to 250 times 
larger cosmological density for this matter, than for ordinary baryonic, what is by 
10 times larger, than it is allowed for dark matter. The contradiction grows to 20 
orders of magnitude, if we take into account the upper bounds on the presence of 
this matter around us. 
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and, being independent on the initial value, [/-hadron abundance decreases 
inversely proportional to time. 

By definition r u = fi/A atom , where A atom is the averaged atomic weight of 
the considered matter and f\ is the initial [/-hadron to baryon ratio. In the 
pregalactic matter this ratio is determined by r 4 from A) Eq.(9), B) Eq.(lO) 
and C) Eq.(12) and is equal to 

10 -10 for the case A, 
3 • 1(T 9 for the case B, (17) 
2 • 10~ 6 for the case C. 



r-4 

n 



Here r b xs 10 10 is baryon to entropy ratio, defined in Eq.(40) in Appendix 3. 

Taking for averaged atomic number density in the Earth n w 10 23 cm -3 , one 
finds that during the age of the Solar system primordial U -hadron abundance 
in the terrestrial matter should have been reduced down to r 10 -28 . One 
should expect similar reduction of [/-hadron concentration in Sun and all the 
other old sufficiently dense astrophysical bodies. Therefore in our own body 
we might contain just one of such heavy hadrons. However, as shown later on 
in Section 4, the persistent pollution from the galactic gas nevertheless may 
increase this relic number density to a stationary much larger (r 10~ 23 ) 
value. 

The principal possibility of strong reduction in dense bodies for primordial 
abundance of exotic charge symmetric particles due to their recombination in 
unstable charmonium like systems was first revealed in [13] for fractionally 
charged colorless composite particles (fractons). 

The [/-hadron abundance in the interstellar gas strongly depends on the mat- 
ter evolution in Galaxy, which is still not known to the extent, we need for 
our discussion. 

Indeed, in the opposite case of low density or of short time interval, when the 
condition (15) is not valid, namely, at 



1 A T tjj _ 3 

n < — ; — r - = A at0 m • 77 • — cm 

r u (av)t 300K t 



4 ■ 10 4 for the case A, 
10 3 for the case B, (18) 
2 for the case C, 



where tjj = A - 10 17 s is the age of the Universe, [/-hadron abundance does not 
change its initial value. 
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In principle, if in the course of evolution matter in the forming Galaxy was 
present during sufficiently long period (t ~ 10 9 yrs) within cold (T ~ 10K) 
clouds with density n ~ 10 3 cm -3 £/-hadron abundance retains its primordial 
value for the cases A and B (f { = f 4 ), but falls down fi = 5- 1CT 9 in the case 
C, making this case close to the case B. The above argument may not imply 
all the ?7-hadrons to be initially present in cold clouds. They can pass through 
cold clouds and decrease their abundance in the case C at the stage of thermal 
instability, when cooling gas clouds, before they become gravitationally bound, 
are bound by the external pressure of the hot gas. Owing to their large inertia 
heavy ?7-hadron atoms from the hot gas can penetrate much deeply inside the 
cloud and can be captured by it much more effectively, than ordinary atoms. 
Such mechanism can provide additional support for reduction of £7-hadron 
abundance in the case C. 

However, according to Appendix 7, in particular annihilation of £/-hadrons 
leads to multiple 7 production. If [7-hadrons with relative abundance f'4 anni- 
hilate at the redshift z, it should leave in the modern Universe a background 
7 flux 

F(E > E 7 ) = N ^f^ r y n ^od-c ^ 3 1Q 3 /4(cm 2 . s . ster) -i ; 

47T 

of 7 quanta with energies E > E 1 = 10GeV/(l + z). Here n imod xs 400cm -3 
is the modern number density of relic photons and the numerical values for 7 
multiplicity iV 7 are given in table 1 of Appendix 7. So annihilation even as early 
as at z ~ 9 leads in the case C to the contribution into diffuse extragalactic 
gamma emission, exceeding the flux, measured by EGRET by three orders of 
magnitude. The latter can be approximated as 

E 

F(E > £ 7 ) « 3 • 10- 6 (-^) L1 (cm 2 • s ■ ster)' 1 , 

Ery 

where E = A51MeV. 

The above upper bound strongly restricts (/ 4 < 10~ 9 the earliest abundance 
because of the consequent impossibility to reduce the primordial [7-hadron 
abundance by £/-hadron annihilation in low density objects. In the cases A 
and B annihilation in such objects should not take place, whereas annihila- 
tion within the dense objects, being opaque for 7 radiation, can avoid this 
constraint due to strong suppression of outgoing 7 flux. However, such con- 
straint nevertheless should arise for the period of dense objects' formation. 
For example, in the course of protostellar collapse hydrodynamical timescale 
t H ~ 1/^/nGp ~ 10 15 s/\/^ exceeds the annihilation timescale 

1 10 12 s /l/30\ 2 T , m 11/10 

m fan{av) fan \ a J l 300iT J [ 250GeV ) 
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at » > 10" (V) 8 m' (sfe)"* (afcr)"*. » * to cm-. We 

consider homogeneous cloud with mass M has radius R ps 1 ° 1/3" (j^) 1 ^ 3 ; 
where M Q = 2-10 33 (? is the Solar mass. It becomes opaque for 7 radiation, when 
this radius exceeds the mean free path Z 7 ~ 10 26 cm/n at n > 3-10 10 cm(-^) 1 / 2 . 
As a result, for / 4 as large as in the case C, rapid annihilation takes place when 
the collapsing matter is transparent for 7 radiation and the EGRET constraint 
can not be avoided. The cases A and B are consistent with this constraint. 

Note that at / 4 < 5- 10~ 6 (-^) 2//9 , i.e. for all the considered cases energy release 
from [7-hadron annihilation does not exceed the gravitational binding energy 
of the collapsing body. Therefore, [7-hadron annihilation can not prevent the 
formation of dense objects but it can provide additional energy source, e.g. 
at early stages of evolution of first stars. Its burning is quite fast (few years) 
and its luminosity may be quite extreme, leading to a short inhibition of star 
formation. 

According to above arguments [7-baryon abundance in the primary cosmic 
rays can be close to the primordial value / 4 . It gives for case B 

/ 4 = r -± rsj 3 • 1(T 9 . (19) 



If (Uuu) is the lightest ?7-baryon, its fraction in cosmic ray helium component 
can reach in this case the value 

( Uuu ) o in -8 

which is accessible for future cosmic rays experiments, such as RIM-PAMELA 
and AMS 02 on International Space Station. 

Similar argument in the case C would give for this fraction ~ 2- 10~ 5 , what may 
be already excluded by the existing data. However, it should be noted that the 
above estimation assumes significant contribution of particles from interstellar 
matter to cosmic rays. If cosmic ray particles are dominantly originated from 
the purely stellar matter, the decrease of ?7-hadron abundance in stars would 
substantially reduce the primary £/-baryon fraction of cosmic rays. 



4 Galactic blowing of £/-baryon atoms polluting our Earth 



Since the condition Eq.(18) is valid for the disc interstellar gas, having the 
number density n g ~ lcm~ 3 one can expect that the Z7-hadron abundance in 
it can decrease relative to the primordial value only due to enrichment of this 
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gas by the matter, which has passed through stars and had the suppressed U- 
hadron abundance according to Eq.(16). Taking the factor of such decrease of 
the order of the ratio of total masses of gas and stars in Galaxy f g ~ 1CT 2 and 
accounting for the acceleration of the interstellar gas by Solar gravitational 
force, so that the infalling gas has velocity v g 4.2 • 10 6 cm/s in vicinity of 
Earth's orbit, one obtains that the flux of £/-hadrons coming with interstellar 
gas should be of the order of 

h = « 1.5 • 10- 7 -^(cm 2 ■ s ■ ster)-\ (20) 



where / 4 is given by the Eq.(17). 

Presence of primordial fJ-hadrons in the Universe should be reflected by their 
existence in Earths atmosphere and ground. However, according to Eq.(16) 
(see discussion in Section 3) primordial terrestrial t/-hadron content should 
strongly decrease due to radiative recombination, so that the ?7-hadron abun- 
dance in Earth is determined by the kinetic equilibrium between the incoming 
L^-hadron flux and the rate of decrease of this abundance by different mecha- 
nisms. 

In the successive analysis we'll concentrate our attention on the case, when the 
lightest U baryon is doubly charged (Uuu) and the lightest U meson is elec- 
trically neutral (Uu). In this case U baryons look like superheavy anomalous 
helium isotopes. 

Searches for anomalous helium were performed in series of experiments based 
on accelerator search [14], spectrometry technique [15] and laser spectroscopy 
[16]. From the experimental point of view an anomalous helium represents a 
favorable case, since it remains in the atmosphere whereas a normal helium is 
severely depleted in the terrestrial environment due to its light mass. 

The best upper limits on the anomalous helium were obtain in [16]. It was 
found by searching for a heavy helium isotope in the Earth's atmosphere that 
in the mass range 5 GeV - 10000 GeV the terrestrial abundance (the ratio 
of anomalous helium number to the total number of atoms in the Earth) of 
anomalous helium is less than 3- 10 -19 - 2- 10~ 19 . The search in the atmosphere 
is reasonable because heavy gases are well mixed up to 80 km and because the 
heavy helium does not sink due to gravity deeply in the Earth and is homo- 
geneously redistributed in the volume of the World Ocean at the timescale of 
10 3 yi\ 
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The kinetic equations, describing evolution of anomalous helium and [/-mesons 
in matter have the form 



-jf = 3u -nu-nfj- {av) - j gU (21) 

for [/-meson number density n and 

= 30 -rifj-nu ■ (av) - j gU (22) 



drift 



dt 



for number density of anomalous helium nu- Here ju and jo take into account 
the income of, correspondingly, [/-baryons and U -mesons to considered region, 
the second terms on the right-hand-side of equations describe UU recombi- 
nation and the terms j g u and 2 g v determine various mechanisms for outgoing 
fluxes, e.g. gravitationally driven sink of particles. The latter effect is much 
stronger for [/-mesons due to much lower mobility of £/-baryon atoms. How- 
ever, long range Coulomb like interaction prevents them from sinking, provided 
that its force exceeds the Earth's gravitational force. 

In order to compare these forces let's consider the World's Ocean as a thin 
shell of thickness L 4 ■ 10 5 cm with homogeneously distributed y charge, 
determined by distribution of [7-baryon atoms with concentration n. The y- 
field outside this shell according to Gauss' law is determined by 

2E y S = Ane y nSL, 

being equal to 

E y = 2ire y nL. 
In the result y force, exerting on [/-mesons 

F = e E 

exceeds gravitational force for [7-baryon atom concentration 

- m /30 _1 \ , , . 



Note that the mobility of ?7-baryon atoms and U mesons differs by 10 order 
of magnitude, what can lead to appearance of excessive y-charges within the 
limits of (23). One can expect that such excessive charges arise due to the 
effective slowing down of [7-baryon atoms in high altitude levels of Earth's 
atmosphere, which are transparent for U mesons, as well as due to the 3 order 
of magnitude decrease of U mesons when they enter the Earth's surface. 
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Under the condition of neutrality, which is strongly protected by Coulomb-like 
y- interaction, all the corresponding parameters for [/-mesons and £/-baryons 
in the Eqs.(21)-(22) are equal, if Eq.(23) is valid. Provided that the timescale 
of mass exchange between the Ocean and atmosphere is much less than the 
timescale of sinking, sink terms can be neglected. 

The stationary solution of Eqs.(21)-(22) gives in this case 

n = fwY (24) 



where 



... 12 fa -I 



ju = 3u=3~^ = W- u - [ ^cm- i s-' (25) 



and (av) is given by the Eq.(3). For j < 10 12 j^ w cm 3 s 1 and (av) given 
by Eq.(3) one obtains in water 




io- 10 

It corresponds to terrestrial £/-baryon abundance 



r< 1Q -23. / fa 



io- 10 ' 



being below the above mentioned experimental upper limits for anomalous 
helium (r < 10~ 19 ) even for the case C with / 4 = 2 • 1CT 6 . In air one has 



For example in a cubic room of 3m size there are nearly 27 thousand heavy 
hadrons. 

Note that putting formally in the Eq.(24) the value of (av) given by the Eq.(6) 
one obtains n < 6-10 2 cm~ 3 and r < 6-10~ 20 , being still below the experimental 
upper limits for anomalous helium abundance. So the qualitative conclusion 
that recombination in dense matter can provide the sufficient decrease of this 
abundance avoiding the contradiction with the experimental constrains could 
be valid even in the absence of gauge y-charge and Coulomb-like y-field in- 
teraction for fJ-hadrons. It looks like the hadronic recombination alone can 
be sufficiently effective in such decrease. However, if we take the value of 
(av) given by the Eq.(7) one obtains n by the factor of ^ ~ 5.5 larger and 
r < 3.3 • 10~ 19 , what exceeds the experimental upper limits for anomalous 
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helium abundance. Moreover, in the absence of ^/-attraction there is no dy- 
namical mechanism, making the number densities of [7-baryons and [/-mesons 
equal each other with high accuracy. So nothing seems to prevent in this case 
selecting and segregating [7-baryons from [/-mesons. Such segregation, being 
highly probable due to the large difference in the mobility of [/-baryon atoms 
and U mesons can lead to uncompensated excess of anomalous helium in the 
Earth, coming into contradiction with the experimental constrains. 

Similar result can be obtained for any planet, having atmosphere and Ocean, 
in which effective mass exchange between atmosphere and Ocean takes place. 
There is no such mass exchange in planets without atmosphere and Ocean 
(e.g. in Moon) and £/-hadron abundance in such planets is determined by the 
interplay of effects of incoming interstellar gas, UU recombination and slow 
sinking of fJ-hadrons to the centers of planets. (See Appendix 6) 



5 Correlation between cosmic ray and large volume underground 
detectors' effects 

Inside large volume underground detectors (as Super Kamiokande) and in 
their vicinity £/-hadron recombination should cause specific events ("spheri- 
cal" energy release with zero total momentum or "wide cone" energy release 
with small total momentum), which could be clearly distinguished from the 
(energy release with high total momentum within "narrow cone") effects of 
common atmospheric neutrino - nucleon-lepton chain (as well as of hypothet- 
ical WIMP annihilation in Sun and Earth). 

The absence of such events inside 22 kilotons of water in Super Kamiokande 
(SK) detector during 5 years of its operation would give the most severe con- 
straint 

n < l(T 3 cm" 3 , 

corresponding to r < 10~ 26 . For the considered type of anomalous helium 
such constraint would be by 7 order of magnitude stronger, than the results of 
present direct searches and 3 orders above our estimation in previous Section. 

However, this constraint assumes that distilled water in SK does still contain 
polluted heavy hadrons (as it may be untrue). Nevertheless even for pure 
water it may not be the case for the detector's container and its vicinity. 
The conservative limit follows from the condition that the rate of [7-hadron 
recombination in the body of detector does not exceed the rate of processes, 
induced by atmospheric muons and neutrinos. The presence of clustered-like 
muons originated on the SK walls would be probably observed. The possibility 
to detect upgoing muon signal from £/-hadron recombination in atmosphere 
is considered in Appendix 8. 
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High sensitivity of large volume detectors to the effects of £/-hadron recombi- 
nation together with the expected increase of volumes of such detectors up to 
lkm 3 offer the possibility of correlated search for cosmic ray fJ-hadrons and 
for effects of their recombination. 

During one year of operation a lkm 3 detector could be sensitive to effects 
of recombination at the [7-hadron number density n ~ 7 ■ lCT 6 cm~ 3 and 
r 7 ■ 10 -29 , covering the whole possible range of these parameters, since 
this level of sensitivity corresponds to the residual concentration of primordial 
[/-hadrons, which can survive inside the Earth. The income of cosmic U- 
hadrons and equilibrium between this income and recombination should lead 
to increase of effect, expected in large volume detectors. 

Even, if the income of anomalous helium with interstellar gas is completely 
suppressed, pollution of Earth by £/-hadrons from primary cosmic rays is pos- 
sible. The minimal effect of pollution by [7-hadron primary cosmic rays flux 
Ijj corresponds to the rate of increase of £7-hadron number density j ~ 
where Re ~ 6 • 10 8 cm is the Earth's radius. If incoming cosmic rays doubly 
charged [/-baryons after their slowing down in matter recombine with elec- 
trons we should take instead of Re the Ocean's thickness L ps 4 • 10 5 cm that 
increases by 3 orders of magnitude the minimal flux and the minimal number 
of events, estimated below. Equilibrium between this income rate and the rate 
of recombination should lead to N ~ jVt events of recombination inside the 
detector with volume V during its operation time t. 

For the minimal flux of cosmic ray fJ-hadrons, accessible to AMS 02 experi- 
ment during 3 years of its operation (I m in ~ 10~ 9 / a ~ 4 ■ 10^ u 1(E), where 
1(E) is given by Eq.(50), in the range of energy per nucleon 1 < E < lOGeV) 
the minimal number of events expected in detector of volume V during time 
t is given by N min ~ ^^" Vt. It gives about 3 events per 10 years in Su- 
perKamiokande (V = 2.2-10 10 cm 3 ) and about 10 4 events in the lkm 3 detector 
during one year of its operation. The noise of this rate is one order and half 
below the expected influence of atmospheric v^. 

The possibility of such correlation facilitates the search for anomalous helium 
in cosmic rays and for the effects of £/-hadron recombination in the large 
volume detectors. 

The previous discussion assumed the lifetime of [/-quarks r exceeding the age 
of the Universe t\j. In the opposite case r <t\j all the primordial U hadrons 
should decay to the present time and the cosmic ray interaction may be the 
only source of cosmic and terrestrial U hadrons (see Appendix 9). 
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6 Signatures for [7-hadrons in accelerator experiments 



The assumed value of £/-quark mass makes the problem of its search at ac- 
celerators similar to the case of t-quark. However, the strategy of such search 
should take into account the principal difference from the case of unstable top 
quark. One should expect that in the considered case a stable hadron should 
be produced. 

If these hadrons are neutral and doubly charged, as we dominantly considered 
above, the last ones can be observed as the 'disagreement' between the track 
curvature (3-momentum) 

p = 0.3B ■ R ■ Q, (26) 



where B is magnetic field in T, p is momentum in GeV, R is radius of curvature 
in meters, Q is the charge of U hadrons in the units of elementary charge e, and 
the energy of the track measured in the calorimeter (or energy loss dE/dx). 
The search for such tracks at the Tevatron is highly desirable. The expected 
inclusive cross section of [/-meson production at the Tevatron is about 0.7 pb 
(for Mq = 250 GeV) and about 0.05 pb for the double charged barion (in the 
case of U-quark); more than a half of the events contains two non-relativistic 
heavy particles with the velocity (3 < 0.7. 

The argument favoring (U uu) as the lightest U baryon is simply based on the 
mass ratio of current u and d quarks. On the other hand there is an argument in 
favor of the fact that the U ud baryon must be lighter than the U uu one. Indeed, 
in all models the scalar-isoscalar -urf-diquark is lighter than the vector-isovector 
uu-diquark. One example is the model with the effective t'Hooft instanton 
induced four quark interaction, which provides a rather strong attraction in 
the scalar w<i-channel and which is absent in the vector uu channel. In a more 
general form we can say that the interaction in isoscalar channel (isoscalar 
potential) must be stronger than that in the isovector case. Otherwise we will 
obtain a negative cross section for one or another reaction since the isovector 
interaction changes the sign under the replacement of rf-quark by the w-quark. 

Thus it is very likely that the Uuu-baxyon will be unstable under the decay 
(Uuu) — > (Uud) + 7r + . This expectation is confirmed by the properties of 
the charmed baryons (where the charm quark is much heavier than two other 
quarks). Indeed, the branching of the S c — > A c + it decay is about 100% [20]. 

The uncertainty in the choice of the lightest baryon of the 4th generation 
is accomplished by the uncertainty of the choice of the lightest quark, is it 
[/-quark in the analogy with the first generation, or D-quark, as it follows 
from the analogy with second and third generations. These ambiguities in 
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theoretical expectations make us to stipulate the signatures for 4th genera- 
tion in accelerator searches, which are independent of the above mentioned 
uncertainties. 

The expected cross section of the charge 1 heavy meson (like Ud) production 
is not too small. It is comparable with the t-quark cross section, depending 
on the mass of heavy quark. At the moment the best limits on the quarks of 
4-th generation was given by the CDF collaboration using the dE/dx mea- 
surements. For the quark with electric charge q = 2/3 the limit is Mjj > 220 
GeV [21]. Therefore in the present paper we choose for our estimates the value 
Mu = 250 GeV, which corresponds to the production cross section of the order 
of 1 pb at the Tevatron energy 4 . 

Due to a very large mass (more than 150 GeV) the new heavy hadrons are 
rather slow. About half of the yield is given by particles with the velocity 
P < 0.7. To identify such hadrons one may study the events with a large 
transverse energy (say, using the trigger - Et > 30 GeV). The signature for 
a new heavy hadrons will be the 'disagreement' between the values of the full 
energy E = ^m 2 + |p| 2 — m measured in the calorimeter, the curvature of 
the track (which, due to a larger momentum \p\ — E//3, will be smaller than 
that for the light hadron where E ~ \p\) and the energy loss dE/dx. For the 
case of heavy hadrons due to a low (5 the energy loss dE ^ D caused by the 
electromagnetic interaction is larger than that for the ultrarelativistic light 
hadron, while in "hadron calorimeter" the energy loss caused by the strong 
interactions is smaller (than for a usual light hadron), due to a lower inelastic 
cross section for a smaller size heavy hadron, like (Ud) meson. Besides this 
the whole large E T will be produced by the single isolated track and not by 
a usual hadronic jet, since the expected energy of the accompanying light 
hadrons E% cc ~ 1 GeV E T is rather low. 

Another possibility to identify the new stable heavy hadrons is to use the 
Cherenkov counter or the time-of-flight information. 

We hope that the hadrons, which contain a heavy quark of 4-th generation, 
may be observed in the new data collected during the Runll at the Tevatron 
and then at the LHC, or the limits on the mass of such a heavy quarks will 
be improved 5 . 



4 The probability to recombine with diquark and to form the baryon is about order 
of magnitude lower. 

5 If the mass of Higgs boson exceeds 2m, its decay channel the pair of stable QQ 
will dominate over the tt, 2W, 2Z and invisible channel to neutrino pair of 4th 
generation [22]. It may be important for the strategy of heavy Higgs searches. 
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7 Discussion 



To conclude, the existence of hidden stable or metastable quark of 4th gen- 
eration can be compatible with the severe experimental constrains on the 
abundance of anomalous isotopes in Earths atmosphere and ground and in 
cosmic rays, even if the lifetime of such quark exceeds the age of the Universe. 
Though the primordial abundance r = r 4 /r fe of hadrons, containing such quark 
(and antiquark) can be hardly less than r ~ 10~ 10 , their primordial content 
can strongly decrease in dense astrophysical objects (in the Earth, in par- 
ticular) owing to the process of recombination, in which hadron, containing 
quark, and hadron, containing antiquark, produce unstable charmonium-like 
quark- antiquark state. 

To make such decrease effective, the equal number density of quark- and 
antiquark-containing hadrons should be preserved. It appeals to a dynami- 
cal mechanism, preventing segregation of quark- and antiquark- containing 
hadrons. Such mechanism, simultaneously providing strict charge symmetry 
of quarks and antiquarks, naturally arises, if the 4th generation possess new 
strictly conserved U(l) gauge (y-) charge. Coulomb-like ^/-charge long range 
force between quarks and antiquarks naturally preserves equal number den- 
sities for corresponding hadrons and dynamically supports the condition of 
^/-charge neutrality. 

It was shown in the present paper that if £/-quark is the lightest quark of the 
4th generation, and the lightest [/-hadrons are doubly charged (?7wu)-baryon 
and electrically neutral (£/-u)-meson, the predicted abundance of anomalous 
helium in Earths atmosphere and ground as well as in cosmic rays is below the 
existing experimental constrains but can be within the reach for the experi- 
mental search in future. The whole cosmic astrophysics and present history of 
these relics are puzzling and surprising, but nearly escaping all present bounds. 

Searches for anomalous isotopes in cosmic rays and at accelerators were per- 
formed during last years. Stable doubly charged U baryons offer challenge for 
cosmic ray and accelerator experimental search as well as for increase of sen- 
sitivity in searches for anomalous helium. In particular, they seem to be of 
evident interest for cosmic ray experiments, such as PAMELA and AMS02. 
+2 charged U baryons represent the low Z/A anomalous helium component 
of cosmic rays, whereas —2 charged U baryons look like anomalous antihelium 
nuclei. However, in the baryon asymmetrical Universe the predicted amount 
of primordial U baryons is exponentially small, whereas their secondary fluxes 
originated from cosmic ray interaction with the galactic matter are predicted 
at the level, few order of magnitude below the expected sensitivity of fu- 
ture cosmic ray experiments. The same is true for cosmic ray +2 charged U 
baryons, if [/-quark lifetime is less than the age of the Universe and primordial 
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U baryons do not survive to the present time. 

The arguments for the lightest (Uuu)-baryon simply use the u and d current 
quark mass ratio. These arguments are not supported by models of quark 
interactions, which favor isoscalar (Uud) baryon to be the lightest among the 
4th generation hadrons (provided that U quark is lighter, than D quark, what 
also may not be the case). If the lightest ?7-hadrons have electric charge +1 
and survive to the present time, their abundance in Earth would exceed the 
experimental constraint on anomalous hydrogen. This may be rather general 
case for the lightest hadrons of the 4th generation. To avoid this problem of 
anomalous hydrogen overproduction the lightest quark of the 4th generation 
should have the lifetime, less than the age of the Universe. 

However short-living are these quarks on the cosmological timescale in a very 
wide range of lifetimes they should behave as stable in accelerator experi- 
ments 6 . 

With all the uncertainties in the predicted mass spectrum of 4th generation 
hadrons one can expect the cross section of their production to be at the 
level of the t-quark cross section. Stability of lightest hadrons and their large 
mass of about 250 GeV should make them rather slow. It offers a chance to 
use as their possible signature the "disagreement" between the total energy, 
measured in calorimeter, on one side, and curvature, single-particle character 
of track and energy loss, on the other side, in the events with large transverse 
energy at Tevatron and LHC. 

In the present work we studied effects of 4th generation having restricted our 
analysis by the processes with 4th generation quarks and antiquarks. However, 
as we have mentioned in the Introduction in the considered approach abso- 
lutely stable neutrino of 4th generation with mass about 50 GeV also bears 
^/-charge. The selfconsistent treatment of the cosmological evolution and as- 
trophysical effects of y-charge plasma of neutrinos, antineutrinos, quarks and 
antiquarks of 4th generation will be the subject of special studies. 

We believe that a tiny trace of heavy hadrons as anomalous helium and stable 
neutral meson 7 may be hidden at a low level in our Universe ~ 10~ 10 — 
10~ 9 ) and even at much lower level here in our terrestrial matter a density 
^ ~ 10~ 23 . There are good reasons to bound the 4th quark mass below TeV 



6 With an intermediate scale of about 10 11 GeV (as in supersymmetry models [25]) 
the expected lifetime of [/-(or D-) quark ~ 10 6 years is much less than the age of 
the Universe but such quark is practically stable in any collider experiments. 

7 Storing these charged and neutral heavy hadrons in the matter might influence 
its e/m properties, leading to the appearance of apparent fractional charge effect in 
solid matter. The present sensitivity for such effect in metals ranges from 10 -22 to 
1(T 20 . 
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energy. Therefore the mass window and relic density is quite narrow and well 
defined, open to a final test. 
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Appendix 1. Freezing out of ?7-quarks 

The expansion rate of the Universe is given by the expression 



V 45 mpi mpi 



which follows from the expression for critical density of the Universe 

_ 3^ _ ^ T 4 
Pcnt ~ 8nG~ 9tot 30 ' 

When it starts to exceed the rate of quark-antiquark annihilation 

Rann = ™4 (<Tv) , (28) 



in the period, corresponding to T = Tj < m, quarks of 4th generation freeze 
out, so that their relative concentration 

r 4 = - (29) 

s 



does not follow the equilibrium distribution at T < Tf. Here the entropy 
density of the Universe on the RD stage is given by 

. = £+£ = |£ = 5^-17- „ 0^.7*. (30) 
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The entropy density can be conveniently expressed in terms of the number 
density n 7 = ^p-T 3 of thermal photons as follows 

S = 45C(3) W7 ~ im 9tot s n T (31) 



Under the condition of entropy conservation in the Universe, the number den- 
sity of the frozen out particles can be simply found for any epoch through the 
corresponding thermal photon number density n 1 . Factors g tot and g tots take 
into account the contribution of all particle species and are defined as 

E( Tbos \ 7 y—- v / Tjer \ 

9bos J + g 2_/ 9fer J 

bosons fermions 

and 

E(Tbos\ 7 ^ (Tfer\ 

9bos \ J + g /_/ 9fer \~^p~ J ' 

bosons fermions 

where gb os and gj er is the number of spin degrees of freedom for bosons and 
fermions, respectively. 

^From the equality of the expressions Eq.(27) and Eq.(28) one gets 

m/T f ^42 + \n(g^ l /2 m p m (av)) 

a 2 

with m p being the proton mass and obtains, taking (av) ~ QC ° and g tot (Tf) = 
gtots(Tf)=g f n 80- 90, 

T f w m/30 

and 

r 4 = — ^— « — ^ ^2.5-10" 14 — . 32 

s f (av) g)l 2 mpi T f (av) 250 GeV 1 ; 



Index "f" means everywhere that the corresponding quantity is taken at T = 
Tf. Also it is worth to emphasize, that given estimation for r 4 relates to only 
4th quark or 4th antiquark abundances, assumed to be equal to each other. 

The modern number density of frozen out particles can be found from this 
estimation as 

^-mod ^4 ' ^mod' 

The modern entropy density s mo d is assumed to be determined by relic photons 
with temperature T and light neutrinos with temperature T v = (4/ll) 1 / 3 T. It 
gives g to t smod = 43/11 and s mod = 7.04 • n imod . 

Note that if Tf > A = m D — m, where m D is the mass of D-quark (assumed 
to be heavier, than [/-quark) the frozen out concentration of 4th generation 



22 



quarks represent at 7/ > T > A a mixture of nearly equal amounts of UU 
and DD pairs. 

At T < A the equilibrium ratio 

D A 

— oc exp (-— ) 

is supported by weak interaction, provided that /^-transitions (U — > D) and 
(£) — > C/) are in equilibrium. The lifetime of D-quarks, r, is also determined 
by the rate of weak (D — > U) transition, and at t 3> r all the frozen out £>D 
pairs should decay to £7?7 pairs. 



Appendix 2. Radiative recombination 



Radiative UU recombination is induced by "Coulomb-like" attraction of U and 
U due to their ^-interaction. It can be described in the analogy to the process 
of free monopole-antimonopole annihilation considered in [9]. Potential energy 
of Coulomb-like interaction between U and U exceeds their thermal energy T 
at the distance 

In the case of y- interaction its running constant a — a y ~ 1/30 [1]. For a < 1, 
on the contrary to the case of monopoles [9] with g 2 /An 3> 1, the mean free 
path of multiple scattering in plasma is given by 

being A ^> d for all T < m. So the diffusion approximation [9] is not valid 
for our case. Therefore radiative capture of free U and U particles should be 
considered. According to [9] , following the classical solution of energy loss due 
to radiation, converting infinite motion to finite, free U and U particles form 
bound systems at the impact parameter 

a w (T/m) 3/10 • d . (33) 



The rate of such binding is then given by 

(av) = 7ta 2 v w 7T • (m/T) 9 / 10 • (-) 2 w (34) 

m 



fi m~ 13 r a \y 300K ^no( 250GeV \nno cm3 

[ l/30 ) T m s 
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The successive evolution of this highly excited atom-like bound system is de- 
termined by the loss of angular momentum owing to y-radiation. The time 
scale for the fall on the center in this bound system, resulting in UU recom- 
bination was estimated according to classical formula in [10] 

a 3 m 2 _ a ( m 21/w 1 



T [ 250GeV } 



As it is easily seen from Eq.(35) this time scale of UU recombination r <C 
m/T 2 <C mpi/T 2 turns to be much less than the cosmo logical time at which 
the bound system was formed. 

Kinetic equation for U-particle abundance with the account of radiative cap- 
ture on RD stage is given by Eq.(4). Its solution at T = T f > T > T QCD = 7\, 
assuming in this period gtot(s) = const — gf what allows to use transformation 
similar to Eq.(41), is given by 

r 4 = PxTiT ^ ^ , ; ~ ~ n^-Z\ T^TTTTTn ( 36 ) 



0.013 ( ™) 1/10 « 4 • IP' 16 ™ ( — 

\Tf) a 2 m P i 250 GeV \ a 



where r is given by Eq.(32). 

At T < Tqcd the solution for the effect of radiative recombination is given by 
r 4 ~ r ° 77777 « r (37) 

with r taken at T = Tq C d equal to r 4 from Eqs.(9),(10) or (12). 

Owing to more rapid cosmological expansion radiative capture of [/-hadrons 
in expanding matter on MD stage is less effective, than on RD stage. So 
the result r 4 « r holds on MD stage with even better precision, than on RD 
stage. Therefore radiative capture does not change the estimation of [/-hadron 
pregalactic abundance, given by Eq.(9). 
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Appendix 3. Hadronization of U quarks 



Here we show that in the baryon excessive background U quarks form U- 
baryons and U form U mesons, while the possible abundance of U antibaryons 
and U mesons is suppressed exponentially. Indeed, even if the number den- 
sity of £7-baryons, n_, was initially comparable with the one of [7-baryons, 
successive reactions with nucleons, such as 

(Uuu) + N -> {Uu) + 2tt (38) 



substantially reduce this number density down to a negligible value. 

The decrease of relative ?7-baryon abundance, r_ = n_/s, is described by the 
equation 

dv 

-7T = - r - ■ n b ■ (<TV) . (39) 



Provided that total baryon number in Universe did not change significantly 
since QCD phase transition, nucleon number density n b at the given time can 
be defined through the relation for baryon to entropy ratio r b , which reads 

W-b , IT-bmod 1 n bmod i n — 10 /,n\ 

— = r b = const = = — — 10 . (40) 

S Smod i -U4 fly mod 

At modern epoch the ratio rjb = nb m od/n Trno d ~ 6 ■ 10~ 10 . 
To solve equation Eq(39) we will use the transformation 

-s-dt = ^^mpi-dT. (41) 



This transformation of differentials is applicable if the effective total number of 
degrees of freedom is approximately a constant during the most of temperature 
interval of considered period following the QCD phase transition, g to t(T < 
Tqcd) = 9tats(T < Tqcd) = Qqcd ~ const. Such a condition can be assumed 
to be fulfilled taking Qqcd ~ 15, which is a roughly averaged value over 
the most temperature interval after QCD transition. This value changes from 
about 17 to about 12 within 150MeV-200keV temperature interval and then 
virtually instantaneously it reduces down to the modern value, which is about 
4. For our estimation we will take the cross section of such recombination as 
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o ~ -\ ■ - and for rate of this reaction (6) 



/ \ 1 a m-^ 011 ^ 

{av) ~ — - ^6-10 



-16 

~ u • ±u 

mf s 



Then one obtains the exponentially strong suppression of primordial [/-baryon 
abundance 



To 

r_ = r exp I — 



■ PHi^H Kro exp(-0.15»^^),(42) 
j s V 45 V m„- mjr/ 



giving for T ~ Tqcd (T\ < T ) 

r_ ~ r exp(-10 10 ). 



Qualitatively the result does not change if we take more conservative estima- 
tion for the rate of this reaction (7). Substituting m p instead of m n in the 
Eq.(42) one obtains for T ~ Tq C d 

r_ ~ r exp (-4 • 10 8 ). 



The same exponential suppression takes place for {Uu) and (Ud) mesons. 
Therefore only [/-baryons and ?7-mesons should be considered as possible pri- 
mordial forms of [/-hadrons. 

Similar suppression of [/-hadrons with light valent antiquarks (u and d) should 
take place in any dense baryonic environment. In particular, it makes the light- 
est U baryon and U meson the only possible forms of ?7-hadrons in terrestrial, 
solar or selenal matter. 

At T > T w pa lMeV the rate of weak /9-reactions exceeds the rate of cosmo- 
logical expansion. Owing to these reactions and small mass difference of light 
u and d quarks one should expect comparable amounts of (Uuu), (Uud) and 
(Udd) baryons (as well as of (Uu) and (Ud) mesons) in the Universe in the 
period Tqcd > T > T w . 

Low number density of [/-hadrons makes negligible their role in reactions 
of Standard Big Bang Nucleosynthesis. In general one may discuss a bound 
deuteron-like state formed by the proton and (Uu) meson in these reactions. 
However it looks unlikely that the interaction between the proton and this 
heavy 'meson' is strong enough to provide such a bound state. 

Indeed, the deuteron is a weakly bound system which has a very small binding 
energy. In conventional proton-neutron potentials (see for example [8]) the 
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attraction between two nucleons is described mainly by the isoscalar exchange 
(such as o"-boson exchange). Therefore, assuming the additive quark model, 
we expect the proton-w-quark interaction to be 3 times smaller than that in 
deuteron. This is not enough to provide the existence of a bound state. 

For the case of p — w-interaction the smallness may be partly compensated 
by the uo exchange which changes the sign interacting with the antiquark. 
Thus the situation with the possibility to form a (Uu) +p bound state is not 
so evident. However, as it was discussed above, after QCD phase transition 
hadronic reactions such as 

(Uu) +p^ (Uud) + 7T° 

should have resulted in exponential suppression of (Uu)-meson number den- 
sity. Moreover, even being formed, (Uu) +p state is unstable relative to this 
reaction. 



Appendix 4. Hadronic recombination 

In this Appendix we evaluate in different ways the cross section for formation 
of charmonium-like (QQ)-meson in collisions of hadrons, containing heavy 
quark Q and its antiquark Q. We assume that such hadrons have initial kinetic 
energy E ~ T, where T is temperature. In the course of collision hadrons form 
"compound hadron" , in which Q and Q- quarks move with relative momentum 
ki n . Heavy quarks also possess ^-attraction with a "fine structure constant" 

Oy. 

In the formation of " compound hadron" y-attraction plays important role. One 
should take into account Sakharov enhancement in the cross section, which at 
v/c<^ 1 is reduced to 2ita y c/v. It reminds the classical problem of accretion, 
when particle falls down the massive body of finite size (radius R) in its cen- 
tral potential. The well known cross section is o ~ 2ttR(R + 2GM/v 2 ). In the 
case of y-attraction the corresponding expression reads a ~ 2nR(R + a y /T), 
with R ~ 1/ m n given by the hadronic size. 

Thanks to the y-attraction two colliding hadrons, one of which contains the 
Q quark and another one - the antiquark Q, obtain the momenta k in ~ 
y/aym^m ~ 1 GeV even for very low T. Thus momentum of heavy quark in- 
side compound hadron is determined by T e ff = max {T, a y m n }. For T < m n 
heavy hadron with momentum /c in crosses the compound hadron (of the size 
R ~ l/m n ) very slowly (with v ~ k in /m) at the timescale t ~ (m/ki n )(l/m w ), 
which is much larger than the hadronization timescale t ~ (l/m n ) So, on one 
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side, relatively large momentum of heavy quarks gives us the possibility to 
neglect the effects of confinement, though, on the other side, the process is 
very slow and hadronization of light hadrons can play important role. 

Since metastable hadrons of 4th generation are not discovered yet (and may 
be do not exist), there is evidently absent any direct experimental information 
about their interactions, and, in particular, about the cross section of reaction 
of hadronic recombination. Moreover, direct information on similar reaction 
for charmed and b-quark hadrons is also absent. It makes this question open 
for theoretical discussions and speculations. Such speculations, assuming for- 
mation of (QQ) -meson as a slow process of successive light hadron evaporation 
lead to the estimation of the order of (6) or (7) for the rate of hadronic recom- 
bination. The argument for such approach is that at low energy of colliding 
hadrons T < m^, it is sufficient to emit few light hadrons (e.g. pions) from the 
compound hadron to bind Q and Q within it and to prevent the disruption of 
compound hadron on (Qqq) and (Qq) states. 

The more justified alternative approach is to take into account relatively large 
momentum of heavy quarks, to neglect effects of confinement and to con- 
sider the forming (QQ) meson as a small size hydrogen-like system. Such 
treatment can provide us realistic minimal estimation of the recombination 
rate. In this approach we can use the well known results for recombina- 
tion of hydrogen atom with the replacement of electron mass m e by the re- 
duced mass M = Mq/2 and the QED coupling a QED by a = Cf^s + oc y ~ 
(4/3) -0.144 + 1/30 = 0.23. 

Here Cf = (N 2 — 1)/2N C = 4/3 is the colour factor; the coupling a s is eval- 
uated at the distances r equal to the size of the ground state. For Mq = 250 
GeV this gives the binding energy E { = (Ma 2 )/(2i 2 ) = 3.2/i 2 GeV 
(Ei = 3.2 GeV for the ground state % — 1. The corresponding momentum 
hi = y/2MEl = 28/i GeV.) 

In this approach the 'inclusive' reaction 

(Qqq) + (Qq) -> (QQ) + ....(anything) (43) 

is considered as the recombination of the "free" heavy quarks 

Q + Q^ (QQ) (44) 

inside the " compound hadron" . It is assumed that since the free quark recombi- 
nation takes place at rather small distances, the influence on it of confinement 
and hadronization of light quarks, taking place at much larger distances with 
the probability equal to 1, can be neglected. 
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Now to evaluate the cross section we can use the known result for the electron- 
proton recombination 

1 87r _ 3 e 4 1 
^W c ¥^ a Mv 2 P (Md 2 /2 + I~) [ ' 



where M and v are the reduced mass and velocity of Q-quark; Ii - ionization 
potential (h = h/i 2 ) [h = E 1 = 3.2 GeV]. The colour factor 1/N C = 1/3 is 
the probability to find an appropriate anticolour. 

To sum approximately over 'i' we note that o r oc 1/i for Ii >> Mv 2 /2 = T e ff 
while at Ii < T e ff the cross section a-i oc 1/i 3 falls down rapidly 



So effectively the sum goes up to i = i max with i max = ^/3.2GeV/T e // corre- 
sponding to I imax = T eff . 

Thus the interpolation formula for recombination cross section reads: 

* = (wr> ' t~T ■ l °z <t^> = L8 ' ' GeV2 (46) 

3 1 T effh T eff T effh 



and the recombination rate is given by 

+ log (T QCD /T eff )) » 0.56 ■ 1<T«> log (Ii/T e ff) ■ ^y/^X^Lyr-^ 

l eff m s 

This formula is valid in the interval of not too small T e ff 3> m 2 r /2M (in order 
to neglect the confinement) and not too large T e ff <C I x (at least T e ff < 

For T e jj = 4.7 MeV (which corresponds to ki n ~ 1.1 GeV) it gives a r = 0.8 
mb and (av) = 7mkbn • c • (250GeV7m) 3 / 2 = 2.1 • io- 19 (252g^)3/2c^ 

Note that here we neglect the momentum of Q-quark inside the (Qqq) or Qq 
hadrons since this momentum k ~ 0.3 — 0.5 GeV is smaller than the incoming 
momentum ki n > 1 GeV, acquired in hadron collision due to y-attraction. 



Appendix 5. Secondary [/-hadrons from cosmic rays 



Astrophysical mechanisms of particle acceleration lead to appearance of charged 
[/-hadron component in cosmic rays. Neutral [/-hadrons can not be acceler- 
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ated in this way directly, but their y-charge can make them to follow accel- 
erated electrically charged [/-particles. Primordial fJ-hadrons can be present 
in the interstellar gas, captured by Solar system. This primary flux of cosmic 
U -hadrons, falling down the Earth should enrich its U -hadron abundance. 

Another source of such enrichment could be cosmic ray interaction with Earth's 
atmosphere. Let's estimate number of quarks of fourth generation that can 
be created in collisions of high energy cosmic protons with nuclei in the at- 
mosphere of Earth. Assume for our estimation that high energy cosmic rays 
dominantly contain protons. To create pair of £/-quark and its antiquark with 
mass m the cm. energy 

y/s > 2m (48) 



is necessary. 

In our case laboratory frame can be connected with Earth. Thus we can esti- 
mate cm. energy of pair proton-nucleon as s ~ 2m p E N , where m p is proton 
mass, E N is the energy of incident nucleon. Then 

2m 2 

E N > E th « w 1.3 • 10 5 GeV. (49) 

m p 



Integral high energy proton flux depends on energy as [12] 



1(E) = 



E 



-1.7, 



cm 2 ■ s ■ ster) 1 , E < E knee , 



k 3 • 10_10 ' (io4ev) ^ ( cm2 ■ s ■ ster)' 1 , E > E knee , 



(50) 



where E knee is about 3 • 10 6 GeV. 

In the case of m — 250 GeV the expected inclusive cross section of £7-hadron 
production in the proton-nucleon collision is of the order of 0.1 pb at E N about 
4.5 • 10 6 GeV and more than 1 pb at E N > 10 7 GeV. To calculate the flux of 
[/-hadrons in cosmic rays we have used the known leading order parton-parton 
cross sections of gg — > QQ and qq — > QQ convoluted with the LO MRST2001 
partons [26] and multiplied the result by the next-to-leading order K factor 
(of about 1.2) and by the mean atomic number of 'Air', A atom = 14.5. The 
ratio of this cross section to the inelastic p — Air cross section, parameterized 
in [23] was convoluted with the proton spectrum Eq.(50) (taken in differential 
form) and integrated over the incoming proton energy. This gives 

Iu = 2.5 -1Q- 21 (cm 2 ■ s- ster)- 1 . (51) 
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About half of the flux ly comes from the region of En > 12.5 • 10 6 GeV where 
the cross section <j\j > 2 pb. The major part of the produced heavy quarks 
after the hadronization forms the Uq and U q mesons. The probability to form 
a baryon in the recombination with a light diquark is expected to be an order 
of magnitude smaller. 

The estimated secondary flux of U -hadrons is by more than 12 orders of mag- 
nitude less, than the primary [/-baryon component of cosmic rays if its fraction 
is given by Eq.(19). 

The production and propagation of heavy hadrons, created by cosmic rays 
at E p > 1.3 • 10 5 ( 25Q ™ eV , ) 2 GeV, will show itself as a long penetrating elec- 
tromagnetic track, whose behaviour will be a combination of lepton-like and 
hadron-like shower. 



Appendix 6. [/-baryon atom pollution in Solar system 



The y-charge neutrality condition holds neutral U mesons near much less mo- 
bile [/-baryon atoms in the slow sinking to the centers of planets. Representing 
the form of anomalous helium atoms [/-baryon atoms are not chemically ac- 
tive and they should not form chemical compositions with matter, which might 
prevent them from sinking down the center. 

Planet's gravity force on the surface is F g = mg and [/-baryon atom of mass m 
sinks due to its action with velocity V = /iF g , where the mobility /i = mnav 
is determined by matter atom number density n and rate of multiple atomic 
collisions with matter av. It gives 

V = —. (52) 
nav 



Taking for Moon g = 1.6-10 2 cm/s 2 ; n = 6-10 23 /Acm 3 ; a ~ 10 16 (m a /m)cm 2 , 
averaged atomic weight in selenal matter A = m a /m p m 30 and thermal 
velocity of atoms v = 2 • 10 4 (m/m a ) 1 / 2 (T/300A') 1 / 2 cm/s, one obtains the 
velocity of sinking down the selenal surface V ~ 10 _8 (300A'/T) 1 / 2 cm/s. 

Inside the homogeneous planet with matter density p at the radius R gravity 
force is F g = mG^-pR. Then 

4tt Gm a R 
3 av 

where m a is the mass of planet matter atom, and the timescale of sinking 
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down the center of planet is given by 



t = R/V = -i ™ = 6 • 10 16 (^— Y' 2 s 

independent of actual matter density and radius of planet. At internal temper- 
ature T > 3000A" this timescale exceeds the age of Solar system, what means 
that sinking can still continue in deep hot planet's interiors. 

Taking for Moon the same flux (20) of [/-hadrons falling down the Moon with 
interstellar gas and taking j ~ ^ E ^ L with R ~ 10 8 cm being the Moon's radius 
one obtains from (24) the number density of [/-hadrons in selenal matter 
n ~ 0.3cm~ 3 . 

In stars [/-baryons should be dominantly ionized and the cross section of their 
collision in plasma does not differ much from such cross section for neutral 
[/-mesons. It strongly increases the velocity of [/-hadron sinking and reduces 
the timescale of sinking down to stellar center. For Sun this timescale is 4 • 
10 9 (T/10 6 iT) 1/2 s. With the account for the order of magnitude increase in 
Eq.(20) of velocity of gas, falling down the Sun, and taking for Sun j ~ 

-TL0 

with Rq = 7 ■ 10 10 'cm one obtains from Eq.(24) the number density of U- 
hadrons inside the Sun n ~ 7 (T '/ 10 6 K ' ' 2 cm~ 3 . Note that the presence of 
doubly charged U-baryons in Solar corona could lead to anomalous He + lines, 
which can be hardly observable if there is no mechanism for local increase of 
U-baryon abundance. 

Effects of Solar wind and Solar activity prevent interstellar gas from falling 
down the Sun reducing Solar abundance of U-baryons. At the distance R from 
Sun such effects decrease as oc R~ 1 / 2 . 

Heliopause at the border of Solar system, where solar wind stops pushing out 
infalling gas, may provide a shield, preventing interstellar gas penetration in 
Solar System. Provided that this shielding is sufficiently effective, anomalous 
helium coming with interstellar gas can be stopped in heliopause, and the 
effects of its pollution of Earth's matter can strongly decrease. 



Appendix 7. Effects of Q-hadron recombination in matter 

Within the effective thickness of atmosphere-Ocean layers or within the sinking 
depth equilibrium between [/-hadron income and recombination is established. 
In the result of recombination, which takes place homogeneously in all the 
regions occupied by U hadrons, rest energy of [/-quarks converts into the 
energy of ordinary hadrons. The latter (dominantly n and K mesons) should 
in their turn give rise to photons and neutrinos as their decay products. Total 
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hadronic energy release in recombination corresponds to 2mlu with energy 
of ordinary hadrons E < m. For Ijj given by Eq.(20) and m = 250GeV it is 
by order of magnitude smaller than the hadronic energy released in Earth's 
atmosphere and ground by cosmic rays. 

Moreover, the dense matter is opaque for charged pions and kaons originated 
from UU recombination. For charged pions with energy E ~ lOGeV the rate 
of their absorption in water 

Y i = n (av) ^4-10 8 l/s 

is by 3 orders of magnitude higher, than their decay rate 

r, = ^-^5-io 5 i/ s . 

E T n 

It leads to strong absorption of charged pions and to corresponding suppres- 
sion (~ Y d /Yi) of neutrino fluxes from their decay. Such fluxes would be still 
larger than the effect of UU recombination in atmosphere, in which it is much 
stronger suppressed by the factor oc (n atm /n w ) 2 ~ 10~ 6 , where n atm and n w 
are respectively atomic number densities in air and water. 

Neutral pions can dominantly escape absorption, but deep layers of water are 
opaque for gamma radiation from their decay, reducing the observable gamma 
source to a surface layer of thickness about Z 7 , mean free path of 7 in water. All 
the Earth is however transparent for neutrino with energy E < m = 250GeV. 

So the stationary regime of UU recombination in Earth should be accompanied 
by the neutrino flux and by the gamma radiation from the Ocean surface. The 
calculation of these fluxes is valid for the general case of any choice for lightest 
Q-quark and Q hadrons, what we reflect in notations below. 

The mean multiplicities of v and 7 produced in the QQ recombination at 
rest were evaluated using the JETSET 7.4 Monte Carlo code [24]. We neglect 
the spin-spin interaction and assume that just due to statistics 25% of the 
UU pairs are in the pseudoscalar (0~) state, and thus mainly decay onto the 
two gluon jets, while 75% of pairs are in the 1~ state and decay onto the 
three gluon jets. The results are presented in Table 1 where we show the total 
multiplicities of gamma and electron, muon and r neutrinos for the case of 
m = 250 GeV. To give some impression about the energy distributions we 
present also the multiplicities of particles with the energy E > 0.1, 1, 10 and 
100 GeV. Besides this in Table 1 we show the energy fractions carried by the 
gamma and electron, muon and r neutrinos as well. Note that the numbers 
for E > 100 GeV had in our calculations a low statistical significance and 
are shown just to demonstrate the strong suppression of the fluxes at high 
energies. Moreover all the numbers corresponds to the QQ annihilation in 
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vacuum and do not account for the attenuation caused by the absorption of 
high energy pions in medium (Earth or Ocean). 

The last effect is crucial for the neutrinos produced via the charged pion decay. 
The lifetime of a fast pion is proportional to the Lorenz gamma factor E/m^ 
and grows with energy. The density of nucleons in the Earth is about n = 
1.5-10 24 cm -3 , the cross section of 7r absorption a(irN) m 20mb = 2-10~ 26 cm 2 
So the mean free path of pions is I = l/(na) = 30cm. On the other hand the 
decay length for charged pions c • r = 780cm • (E/m n ). Thus the suppression 
is l/( CT ) = (30 • U0MeV/780)/E 7T = 5AMeV/E w . Strictly speaking E v < E w 
(so we can expect a stronger than 5MeV/E u suppression) On the other hand 
after the interaction the fast pion does not disappear completely; it creates 
few pions of a lower energy (due to this reason the suppression should be not 
so strong) These two factors approximately compensate each other. So the 
final estimate 5MeV/E u looks reasonable. 

Therefore in the Earth the yield of fast neutrinos will be additionally sup- 
pressed by the factor of about 5MeV/E u . However the prompt neutrinos, 
coming mainly from the decay of charmed particles, are not suppressed by this 
effect. To demonstrate the role of prompt neutrinos for each sort of neutrinos 
in the right columns of Table 1 we present the multiplicities (and correspond- 
ing energy fractions) calculated in the limit when the charged pions, kaons 
and /i-mesons are stable. For a high energies exceeding few GeV the flux of 
prompt neutrinos give the dominant contribution. 





7 


u e + v e 


Vp, + &7i 


v T + v T 


Ntotal 


69 


62 


0.095 


124 


0.091 


0.02 


0.019 


N(E > 0.1 GeV) 


62 


47 


0.094 


97 


0.090 


0.02 


0.019 


N(E > I GeV) 


28 


14.4 


0.080 


31 


0.077 


0.017 


0.016 


N(E > 10 GeV) 


2.4 


0.61 


0.028 


1.5 


0.027 


0.0057 


0.0057 


N(E > 100 GeV) 


0.001 


0.0004 


0.00017 


0.0004 


0.00010 


0.0003 


0.000 


Energy fraction 


0.27 


0.12 


0.0018 


0.26 


0.0017 


0.0004 


0.00035 



Table 1 



Multiplicities of 7, electron-, muon- and r- neutrinos produced in the recombination 
of (QQ) pair with m = 250 GeV. First line shows the total multiplicities while 
the next 4 lines correspond to multiplicities of particles with the energies E > 
0.1, 1, 10, 100 GeV (in the (QQ) pair rest frame). The last line presents the energy 
fraction carried by each sort. For each sort of neutrinos numbers in the left and 
right columns correspond to all neutrinos (produced in the result of fragmentation 
and decays of all annihilating products in vacuum) and to only prompt neutrinos 
respectively. 

In particular, the stationary regime of QQ recombination in Earth should 
be accompanied by gamma radiation with the flux F(E) = N(E)I u l 1 / L, 
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where energy dependence of multiplicity N(E) of 7 quanta with energy E is 
demonstrated in table 1, and Z 7 is the mean free path of such 7 quanta. At E > 
WOMeV one obtains the flux F(E > WOMeV) ~ 3 -lO' 9 ^^ (cm 2 -s-ster)-\ 
coming from the surface layer / 7 ~ 10 2 cm. 

One should also expect neutrino flux from QQ recombination in Earth with 
spectrum given in table. At E v > lGeV it corresponds to the flux ~ 2 • 
10~ 8 T ^ w (cm 2 • s ■ ster)' 1 . 



Appendix 8. Upgoing muon signal from UU annihilation in atmo- 
sphere 

Here we consider the stationary regime of equilibrium between the incoming 
interstellar gas pollution and UU annihilation in matter and discuss the possi- 
bility to detect the upgoing muon signal from such annihilation in atmosphere. 

Assuming an infalling flux, given by Eq.(20) ly ~ 1.5-10~ 7 lc j^ 10 (cm 2 -s- ster)- 1 , 
we must observe a tiny component of annihilation in air (~ 1CT 3 as the ratio 
of atmospheric column height to Ocean's depth). Effects of this annihilation 
will be mostly masked by a more abundant downward cosmic ray showering. 

Netheretheless the presence of a tiny upgoing muons, tracking from UU — > 
hadrons — > — > /i ± chains, is the source of an "anomalous" upgoing muon 
flux. Its value may exceed 

J„ T = N(E fl > 1 GeV)Iu = 4.6 • \^-A_( C m 2 ■ s ■ ster)- 1 . 

This flux is not easy to observe on the ground, but it may be observable from 
mountain, airplane or baloons. 

Indeed, /i trace in air exceeds 6 km its flux estimated above is already com- 
parable with the albedo muon flux, observed at 93°-94° in NEVOD-DECOR 
experiment [27]. Its detection may be easily tested. It may be exceeding 3-4 
orders of magnitude the upgoing muons from atmospheric (7^ ~ (2 — 3) • 
10' 13 (cm 2 ■ s ■ ster)- 1 ). 

Even E^ > 10 GeV harder spectra may play an important role. The expected 
flux is I {E^ > 10 GeV) > 2-10~ lc (cm 2 -s- ster)' 1 with negligible pollution 
by downward cosmic ray secondaries bent by geomagnetic fields. 

These signals could be complementary or even dominant over the expected 
upgoing one from r airshowers due to UHEz/, skimming the Earth's crust 
[28,29,30,31]. 
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Appendix 9. Metastable Q-hadrons 

Interaction of cosmic rays with the spectrum (50) with the matter in disc with 
averaged number density n ~ n g ~ lcm' 3 for the inclusive cross section of 
[7-baryon production taken as in Appendix 5 results during the lifetime of 
galactic cosmic rays T cr ~ 10 7 yr in the creation of integral cosmic ray flux of 
U baryons 

I v w ncr Ub cT cr I (E > Et h ) « 4 • 10" 22 (cm 2 • s ■ ster)' 1 , (53) 

being one order of magnitude smaller than the flux (51) of [/-hadrons produced 
in the cosmic ray interaction with Earth's atmosphere. This flux is several 
orders of magnitude below the level of sensitivity expected in future cosmic 
ray experiments, in particular in AMS 02. For such low incoming flux as (51) 
the number density of anomalous helium stored in Earth during its age is 
n ~ 4 • 10 _12 cm~ 3 , which corresponds to such a low recombination rate that 
kinetic equilibrium is not established and equilibrium solution (24) is not valid. 
The anomalous helium abundance in this case does not exceed r ~ 4 • 10~ 35 , 
what is 16 orders of magnitude below the existing upper limits. It makes in the 
case t < tu the search for anomalous helium at accelerators most promising. 

The above estimations of primordial abundance of [/-baryons and of their 
production rate by cosmic rays are also valid for the case, when the lightest 
[7-baryon is isosinglet (C/i«f) -state with the electric charge +1. However, in 
this case the conclusions differ drastically from the case of double charged 
baryon, since (£/-uci)-baryon represents the form of anomalous hydrogen. Mass 
spectrometry was used on a variety of light elements [14], [17], [18], [19], and 
achieved the lowest limit ~ 10~ 30 on hydrogen isotopic abundance at around 
100 GeV. These experimental constrains are by more than 11 orders of mag- 
nitude stronger, than for anomalous helium. These constrains exclude the 
case of such [/-baryon with r > tu, but in the case r < tu the experimen- 
tal upper limits can be satisfied, since the abundance of anomalous hydrogen 
produced interaction of cosmic rays in atmosphere and stored in Earth is in 
this case about r ps 4 ■ 10 -35 ^-. For a very wide range of r < tu this life- 
time is sufficiently large to prevent the [/-baryon decay in detectors, making 
(meta-)stable [/-baryons a profound signature for 4th generation quarks in 
accelerator searches. 

If particles are sufficiently longliving to survive to the period of galaxy forma- 
tion tg a i, their annihilation can take place and the EGRET constraint on this 
annihilation is valid. For r < t ga i such annihilation can not take place and 
the astrophysical effects of decay should be considered. These effects depend 
on the decay modes, but even for the significant fraction of electromagnetic 
and hadronic cascades among the decay products, effects of [7-hadron decay 
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in CMB spectrum, electromagnetic backgrounds and light element abundance 
are consistent with observational constrains (see for review [32] and [33]) for 
the cases A and B and for a wide range of lifetimes in the case C. In particu- 
lar, for maximal fraction of electromagnetic energy release (~ 0.5) one obtains 
from [33] the constraint f± < 5 • lCr 9 (l + z d ), so that even case C is possible 
for the decay redshift z d > 400, corresponding to the lifetime r < 5 • 10 13 s. 

Similar arguments are valid for the case, when .D-quark is lighter than U- 
quark, and +1 charged (Duu) and (Du) are the lightest D-hadrons. Such 
hadrons can hardly form nuclear bound states with Z > 1 and thus cause the 
problem of anomalous hydrogen overproduction, if .D-quark lifetime exceeds 
the age of the Universe. 

On the other hand, though the current quark mass relation is m d > m u , it is 
not excluded that the self-energy induced by the QED interaction can com- 
pensate the mass difference (m d — m u ), making the electrically neutral (Dd) 
the lightest Z)-hadron. If in the same time the lightest D-baryon is also the 
neutral (Dud), and the lightest D-hadrons are not bound with hydrogen nu- 
clei in anomalous isotopes of hydrogen, primordial D-hadrons would represent 
an interesting form of strongly interacting dark matter particles. Having only 
nuclear interaction with the matter, such particles decouple from plasma at 
the temperatures below T d ~ 30 keV, and then on MD stage, participating 
the development of gravitational instability are dominantly distributed in the 
halo of galaxies. 

After decoupling particle velocity decreases inversely proportionally to the 
scale factor, i.e. as oc T. In this case the rate of radiative binding has the form 

(av) = na 2 v » vr • (m/T d f 10 • (-) 2 (^) 9 / 5 (54) 



and at Trd — T < T < T d — Ti solution (36) of the kinetic equation (4) is 
given by 



r 4 



r 



l + r 



12 



m 



(55) 



\Trd J 



Here Qf = gtotsmod — 43/11 and r is given by (9) in case A, (10) in case 
B or (12) in the case C. In the cases A and B the primordial abundance 
practically does not change, while in the case C the Q-hadron abundance is 
reduced to the end of RD stage by a factor of 10. The energy release, related 
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with the reduction of Q-hadron abundance in the case C, is consistent with 
the constrains on the distortions of CMB spectrum, but hadronic cascades 
from annihilation products can interact with primordial 4 He and pollute the 
primordial light element composition by excessive 6 Li, 7 Li and 7 Be [34] (see 
review in [33]). 



Indeed, following [33], hadronic jet from gluon with energy E g , originated 
from QQ annihilation, contains N p < < ^(Eg/lGeV) 1 ^ 2 pairs of p and 
p with averaged energy E p ~ (Eg/lGeV) 1 ^ 2 . The chain of nuclear reactions 
with energetic 3 He, T, D products of p{p) induced destruction of primordial 
He yields [34] 



4 



(E p /lGeV) 2/3 N p ~ 30 • 2.5 • 10~ 6 U, 



which should not exceed the observed 6 Li abundance 1 ^ M - ~ 10 10 . It restricts 

n b 

the relative amount / 4 of annihilated QQ as / 4 < 1.3 • 10~ 6 , what excludes 
the possibility of QQ annihilation on late RD stage in the case C (for which 
/ 4 = 2-10~ 6 ). Detailed analysis of nucleosynthesis by hadronic cascade, similar 
to the considered recently in [35] would only strengthen this restriction. 

On the MD stage at T RD > T before the period of galaxy formation the rate 
of radiative binding in nearly homogeneously expanding matter retains the 
form (54), but the transformation (41) has now the form 




w/f= ^ ^p,» ( , M ^--dT. (5(3) 



It leads to the solution of the kinetic equation (4) at Trb = T > T > T ga \ = T\ 
given by 



r 4 



1 + r oy I £frnp l J T ° (<rv) v TuD 
r 



(57) 



So on MD stage to the period of galaxy formation at T gal = T CMB • (1 + z gal ) ~ 
20T CMB for T RD = T CMB ■ (1 + z RD ) ~ 2 • 10 4 T CMB (where T CMB = 2.7K is 
the modern temperature of CMB) Q-hadron abundance in the case C should 
decrease by additional factor of 20, making its pregalactic abundance about 
r 4 ~ 10~ 8 . Such Q-hadron annihilation on the MD stage not only leads to 
energy release close to the constrains on CMB spectrum distortions, but also 
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results in the 7 background incompatible with the EGRET data. Thus in the 
case C the considered neutral Q-hadrons should be unstable with the lifetime, 
less than 10 15 s. Note that for the cases A and B, which are not excluded 
by the above arguments, owing to high peculiar velocity effects of Q-hadron 
annihilation in Galaxy are strongly reduced. 

Finally, an interesting case of the lightest (Ddd)-ba.ryon with electric —1 charge 
can be mentioned. At the temperatures about tens keV such baryons can bind 
with protons and helium nuclei in atom-like systems. The ((Ddd) +p)-system 
with the size of ~ 3 ■ 10~ 12 cm can hardly have atomic cross section for its 
interaction with matter. As for the ((Ddd)+He)-system its recombination with 
electrons makes it to follow atomic matter in the period of galaxy formation. 
Specific properties of these systems may provide them to be elusive for the 
existing methods of anomalous hydrogen searches. However, this case is not 
supported by the mass ratios of current quarks or by the arguments of quark 
model. 

Following the additive quark model arguments one can not expect the deuterium- 
like state to be formed by (Dd) with protons. The potential of (Dd)-nucleon 
interaction is expected to be about 3 times smaller than that for p — n inter- 
action. 

For the case of (Dud)-nuc\eon interaction the potential is about 2/3 of the 
p — n potential and the factor (2/3) is not enough to make definite conclusion. 
The problem is that the effective reduced mass M e ff is in this case twice larger 
than that for p — n system, while the true parameter is not the potential V(r) 
but the magnitude V(r) ■ M e ff. So, following the arguments of the Additive 
Quark Model neutral (Dud)n or +1 charged "nuclear" states (Dud)p are not 
excluded and most probably should exist. Their existence would exclude D- 
quark with lifetime exceeding the age of the Universe, even if in the course 
Standard BBN reactions they dominantly transform into Z > 1 states, since 
such transformation can hardly reduce the abundance of Z = +1 states by 15 
orders of magnitude, what is necessary to escape the problem of anomalous 
hydrogen overproduction. So D-quarks with lifetime less than the age of the 
Universe, being hardly eligible to searches for anomalous isotopes in the matter 
and in cosmic rays, can become a challenge for accelerator search. 
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